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Abstract 

In this work, we report a theoretical coupled study of the structural and phonons properties of 
bundled single- and double- walled carbon nanotubes (DWNTs), under hydrostatic compression. 
Our results confirm drastic changes in volume of SWNTs in high-pressure regime as assigned by a 
phase transition from circular to collapsed phase which are strictly dependent on the tube diameter. 
For the DWNTs, those results show first a transformation to a polygonized shape of the outer tube 
and subsequently the simultaneous collapse of the outter and inner tube, at the onset of the inner 
tube ovalization. Before the DWNT collapse, phonon calculations reproduce the experimentally 
observed screening effect on the inner tube pressure induced blue shift both for RBM and tangential 
G z modes . Furthermore, the collapse of CNTs bundles induces a sudden redshift of tangential 
component in agreement with experimental studies. The G 2 band analysis of the SWNT collapsed 
tubes shows that the flattened regions of the tubes are at the origin of their G-band signal. This 
explains the observed graphite type pressure evolution of the G band in the collapsed phase and 
provides in addition a mean for the identification of collapsed tubes. 
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I. INTRODUCTION 



Carbon nanotubes are fundamental materials for developing nanoeletromechanical sys- 
tems because they present a unique combination of outstading electronic and mechanical 
properties. Cross-sectional deformations of single-wall carbon nanotubes (SWNTs) under 
pressure have been extensively studied by means of experimental toolsi"- and theoretical 
models.- - — On the other hand, the behavior of double-wall carbon nanotubes (DWNTs) 
under pressure has been considerably less studied^ - — . DWNTs may be better candidate 
materials than SWNTs for the engineering of nanotube-based composite materials because 
of their geometry, in which the outer tube ensures the chemical coupling with the matrix 
and the inner tube acts as mechanical support for the whole system^. 

Some theoretical studies have been performed to study the pressure dependence of the 
cross-sectional shape of SWNT, depending on the chirality and diameter. By starting from 
an almost perfect circular cross section, the cross- sectional shape of the nanotube becomes 
oval or polygonized as the pressure is increased, evolving later to a collapsed state with a 
peanut shape.— Based on elasticity theory, Sluiter et al. have proposed a diameter- 
dependent phase diagram for the cross-sectional shape of SWNTs under pressure.— Some 
authors reported that the radial deformation of SWNTs for diameters smaller than 2.5nm 
is reversible from the collapsed state while the deformation of larger diameter tubes could 
be irreversible and the collapsed state is metastable or even absolutely stable without pres- 
sure application.— 1 ^ Many experimental studies have suggested that phase transition in 
nanotubes could be dependent on their metallic character or on the surrounding chemi- 
cal environment used for transmitting the pressure.-^— However, theoretical calculations 
suggest that phase transitions of SWNTs under pressure is mainly dependent on the cube 
inverted diameter (p c ~ dj~ 3 ) of the tubes and not on the chirality.— 1 ^ Even if there is a huge 
dispersion of results concerning the pressure transition values from theoretical models, there 
is an overall agreement between different calculations on the existence of two phase transi- 
tions (circular-oval and oval-peanut) and that the critical pressures for those transitions are 
diameter dependent. 

The pressure evolution of DWNTs has been much less studied and, in opposition to 
SWNTs, a detailed understanding of their cross-sectional evolution is still under discussion 
due their complexity regarding the role of inner and outer tubes. Ye et al. has suggested that 
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the critical collapse pressure of DWNT is essentially determined by inner tube stability and 
so the collapse pressure of DWNT is close to what is expected for inner tube^. Other authors 
suggested that collapse pressure of DWNTs is completely different from which is expected 
for the correspondent SWNT when are considerate separately and, the pressure value still 
depends on 1/dJ" 3 scale law but with a suitable choice of an average diameter djf— >2£. In this 
paper we report a study of the vibrational properties of carbon nanotubes bundle under 
pressure, a subject that has not been theoretically well explored in the literature up to now. 

This paper is organized as follows. First, we will describe in section II the methodology 
used to model nanotube bundled structure. In section III. A, we calculated the structural 
evolution of SWNTs and DWNTs bundles under pressure. We focused in structural stability 
of bundle structure by calculating the critical pressure for collapsing and we compare the 
behavior of DWNT bundle with its corresponding SWNT. In section III.B, we explored the 
vibrational properties of SWNT and DWNTs bundle under pressure by calculating RBM 
and tangential phonons modes before and after the nanotube collapse. We close our paper 
with conclusions in section IV. 

II. METHODOLOGY 

In order to access mechanical and vibrational properties of carbon nanotubes under pres- 
sure we initially perform zero-temperature structural minimizations of SWNTs and DWNTs 
bundles. The carbon-carbon bonding within each CNTs is modeled by a reactive empirical 
bond order (REBO) potential proposed by Brenner— 1 ^. Pairwise Lennard- Jones potentials 
are added to model the non-bonding van der Waals terms (e/fc&=44K, a = 3.39 A), which 
are essential to describe intertube interactions^. We studied zigzag and armchair C bundles 
of CNTs in triangular arrangements within orthogonal unit cells containing two SWNTs or 
DWNTs each, with lateral lattice constants a x and a y = y3a x (the a y /a x ratio is kept fixed). 
The lattice constant along the z axis, (a z ), is chosen to contain 5 (8) unit cells of the zigzag 
(armchair) tubes. 

We perform a sequence of small and controllable steps of unit cell volume reduction for 
each fixed nanotube phase studied (circular, polygonized, oval, peanut, etc.). Each step is of 
the order of AV/Vq = -0.01% , where Vq is the ambient pressure volume for each nanotube 
system in circular phase. For each fixed volume V, we search for the atomic positions and 
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lattice constants a x and a z that minimize the internal energy U(V). The pressure is obtained 
by p = —AU / AV and the enthalpy is H = U + pV. 

Phonon frequencies and eigenvectors are directly obtained by the diagonalization of the 
force constant matrix. The matrix elements are calculated by using finite difference tech- 
niques. We focus our vibrational analysis on the radial breathing mode (RBM) and longitu- 
dinal G-band (denoted G z ), whose displacements are along the tube axis. Therefore, RBM 
and G z projected density of states (PDOS) are constructed by projecting the phonon eigen- 
vectors onto the corresponding radial and axial (with opposite phases for atoms in different 
sub-lattices) displacement fields, respectively. 

III. SWNTS AND DWNTS UNDER PRESSURE: STRUCTURAL AND VIBRA- 
TIONAL PROPERTIES 

A. Nanotube Collapse 

We start by studying the low-pressure behavior of SWNT bundles, considering four basic 
cross-sectional shapes namely: circular, oval, polygonal (hexagonal) and peanut-like. Tubes 
of different diameters exhibit a different sequence of cross-sectional shapes as a function of 
pressure^. In particular, the polygonal or hexagonal phase is typical of bundled tubes and 
it is often obtained in calculations for large diameter nanotubes since plane-parallel facing 
between adjacent tubes tends to decrease the van der Walls interaction energy^"—. Fig[TJi 
shows the enthalpy as a function of pressure calculated for circular, oval and polygonized 
phases for bundles of (8,8), (18,0) and (24,0) SWNTs. The crossing of two curves in this 
plot by following the lowest enthalpy for a given pressure indicates a transition between 
two different cross-sectional shapes. The derivative discontinuities in the H(p) plots are 
associated with the volume variation at the phase transition. This can be seen more clearly 
when p — V plots are constructed (not shown here) where each transition is marked by a 
discontinuous change in volume. For the (8,8) SWNT bundle, the transition pressure from 
circular to oval (so called pi) occurs close to 1.5 GPa and from oval to peanut (so called 
p 2 ) around 2.6 GPa. For the (18,0) SWNT bundle, we found phase transitions at 1.2 GPa 
(circular to polygonal, p' x ) and 1.5 GPa (polygonal to peanut, p' 2 ). Finally, for the (24,0) 
SWNT bundle the polygonal to peanut transition (p^) is found at 0.6 GPa. We see that 
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the small diameter (8,8) SWNT bundle undergoes a circular— )-oval—)-peanut phase transition 
sequence as the pressure increases, whereas for the intermediate diameter (18,0) SWNT the 
sequence is circular— j-polygonal— ?-peanut and for the large-diameter (24,0) we obtain simply 
a polygonal-peanut transition, since the tubes are already polygonized even at zero pressure. 




1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 

Pressure (GPa) Pressure (GPa) Pressure (GPa) 



FIG. 1. Enthalpy vs pressure curvers calculated for (a) (8,8), (b) (18,0) and (c) (24,0) bundled 
SWNTs. The pressure evolution for (10,0) SWNT is similar to (8,8) SWNT showed in (a) and the 
circular-oval and oval-peanut transitions is found close to 1.55 and 9.6 GPa 

As expected, the critical pressures are strongly diameter-dependent^^ 1 ^ 1 ^ 1 ^ 1 ^. Actu- 
ally, many authors have pointed out that the critical pressure (pi) for circular-oval transition 
scales with dj~ 3 , where d t is tube diameter.— Some authors also found that the critical pres- 
sure for oval-peanut transition (p 2 ) also scales with d^ 3 .— Our results for bundled SWNTs 
shown in Fig [TJ point out that the transition to the peanut geometry (p 2 and p 2 ) follows ap- 
proximately the same law p 2 = Cjd\ with C=4.4nm 3 .GPa. This scaling law was observed 
for all studied zigzag and armchair bundled tubes. However, the pressure transition values 
Pi e p\ for circular to intermediate phase (oval and polygonal, respectively) show a slightly 
dependence on the inverted tube diameter but they do not follow the djr 3 scaling law. 
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FIG. 2. p-V curves (a) and enthalpy functions (b) for the pressure-induced sequence of phase 
transitions in the (10,0)@(18,0) DWNT. Different cross-sectional shapes correspond to different 
colors in the plots and are shown in the insets. 

We also have modeled several DWNT bundles under pressure and choose the (10,0)@(18,0) 
DWNT (This notation means that the (10,0) tube is inside the (18,0) tube) in order to com- 
pare its structural stability with those of their constituent SWNTs. We find a very similar 
behavior for the other DWNTs such (12,0)@(20,0) and (11,0)@(19,0). In Fig. |2Ji and[2b we 
show respectively the p — V and enthalpy curves for the (10,0)@(18,0) DWNT bundle. Upon 
compressing the original circular structure of a DWNT we find four distinct structures as 
shown in Fig. [2^. The A configuration is the non-deformed structure with a circular shape 
for both inner and outer tube. The B configuration is characterized by the polygonization of 
the outer tube while the inner tube maintain its circular cross section. The C configuration 
corresponds to a polygonal outer tube and oval inner tube. Finally, the D configuration 
stands for the case where the outer and inner tubes display oval/peanut cross- sectional 
shapes. 

The polygonization of outer tube (A- B phase transition) in (10,0)@(18,0) DWNT bundle 
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is observed at 0.81 GPa which is lower than the value (1.2 GPa) what was calculated for 
the same transition in the (18,0) SWNT bundle. This result suggests that the presence of 
the inner tube enhances the outer-outer tube interactions within the bundle, thus reducing 
the critical pressure value for polygonization. Upon increasing the pressure, a B-C transfor- 
mation occurs at 6.2 ± 0.2 GPa, since the inner tube ovalizes and the outer tube reaches a 
higher degree of polygonization. We are able to increase the pressure in the C phase up to 
13 GPa without inducing any strong deformation. However, from the enthalpy analysis we 
conclude that the D phase is more stable than C phase for the whole range of investigated 
pressures. Thus the C configuration is actually metastable. By further increasing the pres- 
sure in the D phase, it reaches the same enthalpy as the B phase around 5.7 ± 0.2 GPa. 
The D phase is also metastable even at lower pressures, down to 5.5 GPa. In summary, the 
actual sequence of phase transitions for the (10,0)@(18,0) DWNT bundle is A — > B at 0.81 
GPa and B — > D at 5.7 ± 0.2 GPa. This value contrasts with the 1.4 GPa calculated col- 
lapse pressure of the (18,0) SWNT outer tube, cleary poiting out to the structural support 
given by the inner tube, as already proposed by some authors^ 1 ^. 

We conclude that the C configuration (polygonized outer tube and ovalized inner tube) 
is unstable when we compare with the both collapsed tubes which let us propose that the 
ovalization of the inner tube is the reason for the overall collapsing of DWNT. For isolated 
DWNTs, Ye et al^ suggest that the DWNT transition could be essentially determined 
by the inner tube transition. They also found that collapse of DWNTs follow almost the 
same dj~ 3 law when the inner tube is used, which means that p c for DWNTs is the same 
as expected for SWNT inner tub©2S. In such way, our results agrees in the fact that the 
inner tube ovalization determines the critical pressure for collapsing of DWNTs. However, 
there is a clear pressure screening effect from the outer tube in the behavior of the inner 
tube because the critical pressure for the expected circular-oval transition for inner tube 
was increased to 5.7 ± 0.2 GPa when it was encapsulated in the DWNT. This value is 
considerably higher than that observed for the (10,0) SWNT bundle, which was 1.55 GPa. 
Second, the full collapsing of the inner tube (oval-peanut transition) is continuous in the 
DWNT case, differently from the SWNT case, in which there is a discontinuous change in 
the volume (around 9.6 GPa). For the DWNT case, this transition is marked by a change 
in compressibility (which is proportional to the slope of the p — V plot) at around 10 GPa, 
clearly seen from Figj2j We found a similar behavior for (12,0)@(20,0) DWNT bundle, where 
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the transition B — > D is predicted to occur at 5.5 ± 0.1 GPa and change in bulk modulus 
takes places around 5.9 GPa close to that value expected for (12,0) SWNT (5.4 GPa). 



B. Phonon Density of States evolution under pressure 
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FIG. 3. Phonon density of states (DOS) projected in RBM (a) and tangential (b) probing vectors 
for (18,0) SWNT. Some snapshots of the bundle structure evolution can be followed in (c). 

Fig |3] shows the phonon projected density of states (PDOS) calculated for (18,0) SWNT 
bundles. The PDOS for RBM and G^ modes are shown in red (a) and blue (b), respectively. 
From the figure, we can follow the RBM and G z evolutions as pressure is increased, and the 
corresponding snapshots of the bundle structural evolution are shown in the right panels. 
First, we observe the RBM mode of (18,0) SWNT in the low frequency region centered at 
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around 150 cm -1 in the circular phase (0.0 GPa), which gradually shifts to higher frequencies 
as the pressure is increased. The pressure coefficient of this mode is 7.0 ± 2.5 cm _1 /GPa, 
in good agreement with experiments^ 1 ^ 1 ^. After collapse, the contributions to the radial 
mode spread in the low frequency region, which is equivalent to say that the RBM is not a 
single, well-defined mode any longer in good correspondance with the experimental difficulty 
to observe RBM after the some suggested pressure phase transitions.-'^L^i-— 

The G z band behavior is shown in the right panel. We clearly observe the splitting of this 
mode when polygonization takes place. After collapse, we clearly see a sudden jump to lower 
frequencies and intensity enhancement of the lower-frequency contribution, while the higher- 
frequency one continues to upshift. We studied several SWNTs bundles and similar results 
were observed. For (10,0) SWNT bundle the spread of radial contribution were observed in 
circular-oval transition i.e. before the collapsing of tube. The split of G z component due 
polygonization is also observed for (24,0) SWNT bundle and, evidently, not observed for 
(10,0) SWNT. Furthermore, the sudden red-shift of G z component after collapsing is also 
observed for (24,0) SWNT. 

The phonon evolution with pressure for (10,0)@(18,0) DWNT bundle can be observed 
in Fig HI It is clear from left panels the presence in low-requency region of two distinct 
peaks around 150 cm -1 and 325 cm -1 corresponding to RBM modes of outer and inner 
tubes, respectively. The peaks are slightly shifted from the corresponding SWNTs in the 
circular conformation, as expected from the intertube coupling^. It is interesting to note 
the evolution of both RBM peaks with pressure, where we clearly observe that the RBM 
frequency of the outer tube shifts much faster than the RBM frequency of the inner tube 
before collapse. This is a clear evidence of a screening effect on the inner tube by outer tube. 
Pressure screening effects on RBM and G band in DWNTs have been observed in several 
Raman experiments^ 1 ^^— . Our calculations of RBM and G z confirm that this pressure 
screening effect occurs well before any structural collapse. In the (b) panels of Fig. HJ we 
can follow the G 2 component for outer and inner tube. It should be noted here that inner 
(10,0) tube G 2 component is well overestimated being at about 1850 cm -1 compared with 
the outer (18,0) one which is located at around 1650cm -1 . However, qualitative analysis 
could be still performed. It is interesting to note that pressure screening effects are also 
observed in tangential components as we can see that the outer tube G 2 components shift 
faster than the inner ones. As observed for SWNTs, the transition of the outer tube to a 
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polygonized phase at 0.85 GPa is clearly marked by a splitting of the G z band. After colapse, 
the PDOS spreads over a large frequency range for radial and tangential contributions but 
it is possible to see that the lowest frequency components of G z are sudden shifted to lower 
frequencies, as in the case of SWNTs bundles. 
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FIG. 4. Phonon density of states (DOS) projected in RBM (a) and tangential (b) probing vectors 
for (10,0)@(18,0) DWNT. Some snapshots of bundle evolution can be followed in (c). 



In Fig. [51 we plot the z-displacement i.e., along the tube axis direction, of some eigen- 
vectors as a function of the coordinate angle 9 defined from the center of (18,0) tube. First, 
we have identified the A\ g mode centered around 1721.4 cm -1 for the structure at -1.0 GPa 
(Fig. [5^). This is the mode which mainly contributes to v-PDOS as we can see in Fig. 
[5b. As the pressure is increased up to 1.4 GPa, the (18,0) cross section is poligonized as we 
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observe a split of G z PDOS contribution. Analysis of the phonon displacements (Fig. [5^) 
shows that the higher frequency component is mainly due to modes which are localized on 
the high-curvature regions (vertices of the hexagons). Modes centered at 1771.3, 1780.7 and 
1786.1 cm -1 have the same symmetry and the their maximums of displacement are localized 
for 9 equal to ±30°, ±90° and ±150° which are the vertices of polygonized shape. After 
the collapse of (18,0) SWNT bundle, we have calculated phonon eigenvectors for collapsed 
structure at 2.2 GPa. In Fig. [5^, we show modes centered at 1800.6, 1803.2 and 1815.4 
cm -1 which mainly contribute to higher-frequency peak in PDOS. Then, we clearly see that 
high-frequency contribution in PDOS at 2.2 GPa come from modes that are localized in the 
high-curved regions of the peanut-shaped tube. {9 ~ ±180° and 9 ~ ±0°). Consequently, 
the low-frequency and more intense peak of PDOS is due to vibration modes which are 
localized in the flat regions of the peanut-shaped structure. 



From an experimental point of view, it has been recently proposed that a saturation 
or even a negative pressure slope of Raman G + component for SWNT and DWNT could 
be associated to the collapse of the nanotubes^^ 1 ^. Futhermore, after the collapse this 
band follows the graphite pressure evolution as observed in Fig. |6j with a smaller pressure 
coefficient^^.. This was also observed in PDOS of SWNT bundle calculations (cf. Fig 
|3b). Our calculations confirm this hypothesis and identifies the atomistic origin of the low- 
frequency bands. Collapsing of the nanotube cross-section leads to flattened regions where 
the reduced stress in C-C tangential bonds reduce the G z frequency. However, as we observe 
in our calculations, a small high-frequency component of G 2 component, arising from the 
high-curvature regions, remains after the collapse. Since the flattened portion of the tubes 
is much higher than curved portion, we expect that the low-frequency shift of the G band 
is dominant in the experiments. The experimentaly red-shift observed of the G-band is 
then explained from our calculations as related to the dominance of the flat region of the 
collapse phase in the Raman signal. This is in addition perfectly coherent with the fact in 
the collapsed state the pressure evolution of the G-band clearly matches (see Fig. [6]) the 
one of graphite or graphene under triaxial compression^^. 
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FIG. 5. (a) Amplitude of z-displacements of (18,0) SWNT bundle as a function of coordinate angle 
before the collapse (-1.0 GPa and 1.4 GPa) and after collapse (2.2 GPa). From bottom to top: 
Eigenvectors for A\ g mode centered at 1721.4 cm^ 1 (black) at -1.0 GPa; for modes centered at 
1771.3, 1780.7 and 1786.1 cm" 1 (red) at 1.4 GPa, and for modes centered at 1800.6, 1803.2 and 
1815.4 cm" 1 (blue) at 2.2 GPa. (b) v-PDOS projected in z direction for (18,0) SWNT at -1.0, 1.4 
and 2.2 GPa. Arrows mark the center of modes whose eigenvectors are plotted in (a). 



IV. CONCLUSIONS 



We studied high-pressure structural modifications and phonon-mode shifts of single and 
double wall carbon nanotubes bundles using zero-temperature enthalpy minimization with 
classical interatomic potentials. We confirmed that the structural transformation of SWNTs 
in bundles from circular to collapsed cross-section under hydrostatic pressure is strongly 
diameter-dependent. It is also evident that the transition from polygonized to collapsed 
cross section for large diameter tubes has a first-order character with a large bundle volume 
discontinuity whereas for small diameters the transition is continuous, with intermediate oval 

12 



1720 



n — i — i — rn — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



1700 



-1680 h 
E 

£ 1660 - 
.c 

c 

E 1640 

or 
+ 

O 1620 



1600 



1580 



l 
» 



J I I l_l I I I I I I I I I I I I I I I I I I I I I I I L 



• SWNT 

• DWNT (outter) 
— Graphite 

L 



_L 







10 15 20 
Pressure (GPa) 



25 



30 



FIG. 6. Experimentaly pressure Raman shift evolution of G + component for SWNTs, DWNT and 
graphite. Adapted from Ref i 1 7 i 39 i 4 'i 48 



or racetrack cross sections. For DWNT bundles, screening pressure effects were observed 
on the inner tube, which keeps its circular cross section for higher pressures than would 
be expected from the SWNTs behavior. Furthermore, the inner tube acts as structural 
support to outer tube. Phonons calculation also reveals screening effects in radial and 
tangential component as the pressure increases. Polygonization of (18,0) SWNT bundle is 
characterized by a split of G z PDOS contribution which is mainly due to localized modes on 
high-curved region of polygonal shape. After the collapse transition, the tangential modes 
associated with the flat regions of the tubes are suddenly shifted to low frequencies and high- 
curved contribution of peanut shape is less pronounced but still observed in G z PDOS. The 
experimentaly observation of the G-band red-shift at nanotube collapse let us conclude on 
the dominance on the Raman signal of collapsed tubes of the flat regions of collapsed tubes. 
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Moreover, this provides a coherent explanation of the lower pressure slope (~ 4cm _1 /GPa) 
of the G-band observed after nanotubes collapse. Consequently, the pressure slope could be 
a useful means for the identification of the geometry shape state of SWNT and DWNTs. 
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